Proximity of antiferromagnetism and superconductivity in LaOi-^F^FeAs: effective 

Hamiltonian from ab initio studies 
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We report density functional theory calculations for the parent compound LaOFeAs of the newly 
discovered 26K Fe-based superconductor LaOi-^F^FeAs. We find that the ground state is an 
ordered antiferromagnet, with staggered moment about 2.3/xs, on the border with the Mott insu- 
lating state. We fit the bands crossing the Fermi surface, derived from Fe and As, to a tight-binding 
Hamiltonian using maximally localized Wannier functions on Fe 3d and As 4p orbitals. The model 
Hamiltonian accurately describes the Fermi surface obtained via first-principles calculations. Due 
to the evident proximity of superconductivity to antiferromagnetism and the Mott transition, we 
suggest that the system may be an analog of the electron doped cuprates, where antiferromagnetism 
and superconductivity coexist. 
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PACS numbers: 74.70.-b,74.25.Ha,74.25.Jb,74.25.Kc 

The recent discovery of superconductivity with onset 
temperature of 26K in LaOi-^F^FeAsi has generated 
considerable interest because of a number of unusual as- 
pects of this material. First, with the exception of some 
of the A-15 materials, Fe is never found in superconduc- 
tors at zero pressure (although Fe itself superconducts at 
lOGPa). Second, both ferromagnetic and antiferromag- 
netic fluctuations are apparently present in the material, 
suggesting possible analogies to ternary rare earth, heavy 
fermion, borocarbide, ruthenate and cuprate supercon- 
ductors. Finally, the discovery of superconductivity at 
the relatively high critical temperature of 26K implies 
that a new pairing mechanism may be in play. The ana- 
log system LaOi-^F^FeP has a critical temperature of 
7K, so there appears to be a new class of superconducting 
materials with no obvious limit on T c . 

Naively, ferromagnetic order is inimical to supercon- 
ductivity since the exchange field of the ferromagnetic 
ion breaks singlet pairs. Ferromagnetic order was there- 
fore found only very recently to coexist with supercon- 
ductivity in UGe^ and URhGa 3 , in situations where the 
order is quite weak. Coexistence of antiferromagnetic 
order, on the other hand, is less pairbreaking if the co- 
herence length is much larger than the wavelength of the 
magnetic modulation, which is generally the case. Thus 
many examples of superconductors coexisting with anti- 
ferromagnetic order are known, and have been recently 
reviewed^. 

The structure of the new material consists of LaO lay- 
ers sandwiching a layer of FeAs, and doping with F ap- 
pears to occur on the O sites. Since the Fe is arranged 
in a simple square lattice, the analogy with the cuprates, 
where electrons hop on a square lattice and doping occurs 
via a nearby oxide charge reservoir layer, is tempting to 
draw. Early electronic structure calculations for both the 
P— and the new As materials&i£ have presented a some- 
what different picture, however. While experimentally a 
low charge density was measured for this material^, these 
calculations suggest that five bands, of primarily Fe-As 



character, cross the Fermi level and give rise to a multi- 
sheeted, quasi-2D material. No evidence for long-range 
order was found in these studies, although proximity to 
both antiferromagnetic and ferromagnetic ordered states 
was noted. 

The weak coupling of the LaO layers to the FeAs lay- 
ers found here and in previous works also suggests that 
insight may be gained by examining iron monoarsenide 
FeAs, a layered metallic helimagnetiS with large pitch 
angle in the zincblende structure with low-temperature 
effective moment ~ 0.5[1b- The electronic structure and 
spin moment for this compound have been calculated by 
density functional theory--. While the theory is success- 
ful in the sense that an antiferromagnetic ground state is 
found, it does not distinguish the complicated magnetic 
structure and finds a moment of order ~ 2^ig. 

To understand the electronic structure properties of 
the new Fe-based superconductors and the interplay 
with structure and magnetic states, we have performed 
first-principles density functional theory simulations on 
the undoped and x=0.0625 as well as £=0.125 doped 
LaOi-^F^FeAs. Most of the reported results on elec- 
tronic structure and model Hamiltonian have been ob- 
tained using the PWSCF packaged, which employs a 
plane- wave basis set and ultrasoft pseudopotentials 13 . 
We have also used VAS P 14 i 15 to confirm our calcula- 
tions when the same calculations can be performed, as 
described in the following sections in detail. The local 
spin density approximation (LSDA) and generalized gra- 
dient approximation (GGA) of Perdew, Burke and Ernz- 
erhof (PBE) 16 potentials have been incorporated for the 
simulation. For density of states (DOS) calculations, we 
have used a 16x16x8 Monkhorst dense gridil to sam- 
ple the Brillouin zone; while for structural relaxation and 
self-consistent calculations, a 8x8x4 Monkhorst grid has 
been used. All structures have been fully optimized un- 
til internal stress and forces on each atom are negligi- 
ble. Our GGA+U calculations have been performed via 
the VASP code. The existing literature uses an on-site 
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Columb energy that varies from 4.0 eV to 6.9 e V 18 i 19 ' 20 , 
and we have explored the parameter space within the 
range U=2. 0-5.0 eV and J=0.89 eV on Fe in our calcu- 
lations. 
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FIG. 1: (Color online)Undoped LaOFeAs band structure 
of (a) AFM state and (b) PM state. Red lines represent 
DFT calculation results; blue lines are band structures re- 
constructed from the tight-binding model using maximally 
localized Wannier functions (MLWF). Since spin- up and spin- 
down bands are degenerate for AFM state, we plot only spin- 
up bands here. For both figures, Fermi energies are indicated 
by the green line at eV. 

Our calculations show an unambiguous antiferromag- 
netic (AFM) ground state with staggered moment 2.3 
Hb for undoped LaOFeAs, which is 84 meV per Fe lower 
than paramagnetic (PM) and ferromagnetic (FM) states. 
The energy difference between the latter two is found to 
be negligible. In fact, the FM state has a very small 
magnetic moment (~0.05 per Fe); therefore it can be re- 
garded as a PM state. The AFM ground state has been 
confirmed by independent VASP calculations using the 
projector augmented wave (PAW) method 2 - 1 . The op- 
timized structure has a lattice constant of a=4.0200 A 
and c=8.7394 A; and the bond lengths for Fe-As and 
La-0 are 2.35 A and 2.40 A respectively. For reference, 
the paramagnetic state has an optimized lattice constant 
of a=3.9899 A and c=8.6119 A, while the bond lengths 
for Fe-As and La-0 are 2.34 A and 2.33 A, respectively. 
Both the AFM and PM band structures are shown in 
Fig. [T] with red curves. In both states, a small disper- 
sion along the c-axis (from T to Z and from A to M) 



indicates interactions between layers are weak, and thus 
the separation of the structure into LaO and FeAs lay- 
ers is possible. The PM state band structure reproduces 
previous DFT calculation results^, exhibiting 5 bands 
across the Fermi level. The AFM state band structure 
is qualitatively different, exhibiting only 3 bands across 
Ep. In VASP calculations, AFM states are 14 meV per 
Fe lower than PM and FM states. Similar to PWSCF cal- 
culations, the band structures of the AFM state are very 
different from the PM state. These results indicate the 
delicacy of the magnetic states in this system, and that 
the magnetism strongly affects the electronic structure. 
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FIG. 2: (Color online)DOS (top panel) and PDOS (FeAs and 
LaO planes: middle panel and bottom panel, respectively) of 
undoped LaOFeAs. Since the spin-up and spin-down states 
are degenerate for AFM state, we plot spin-up states only. 
The Fermi level is aligned to 0.0 eV. 

To further examine and confirm these findings, we have 
performed two series of additional calculations. First, we 
have performed GGA+U calculations using VASP. Haule 
et. al used a dynamical mean field theory (DMFT)-LDA 
approach, and found that a critical value of £7=4.5 eV 
led to a Mott transition with a gap at the Fermi surface. 
We have calculated the electronic structure using VASP's 
implementation of GGA+U, and also find a Mott tran- 
sition for the LaOFeAs system at a critical U ~ 3 eV 
for Fe. Note that lower bound of the empirical value of 
U chosen in calculation is 3.5-4.0 eV for Fe d orbitals 
(_LS). rpj^ g rounc j s t a t e is found to be always AFM for 
all tested U within 0.0-5.0 eV in our calculations, but 
the DOS changes dramatically (Fig. [4|. A Mott gap 
of about 1.0 eV is observed in the GGA+U calculation 
at U — 4.5 eV. Experimentally, it is observed that below 
100K, the resistivity of undoped LaOFeAs increases when 
temperature decreases, but appears to remain metallic 1 , 
suggesting that the system is in fact on the edge of a 
Mott transition.—. 

Second, we have investigated bulk FeAs. The AFM 
state is again found to be the ground state, in agreement 
with experiment^, as well as with previous calculations 
based on full-potential linearized augmented planewave 
(FLAPW) method 11 . In addition, an isolated layer of 
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FeAs has also been simulated, and the system again has 
an AFM ground state. The calculated AFM states of 
bulk and the isolated layer are 53 meV and 408 meV per 
Fe atom lower than their PM states, respectively. Com- 
pared to FLAPW calculations, our results have shown a 
smaller energy difference between AFM and PM state, 
indicating that our calculations do not have an artificial 
bias for the AFM state. 

The electron DOS derived from the AFM band struc- 
ture and the corresponding projected DOS (PDOS) onto 
LaO and FeAs planes are presented in Fig. [5] It is clear 
from the PDOS analysis that Fe 3d orbitals dominate the 
DOS around Ep (-2 to 2 eV relative to Ep) and from -12 
to -10 eV; whereas DOS below -13 eV is almost com- 
pletely derived from LaO layers. This clear separation 
in the DOS confirms that this material can be viewed 
as a layered structure. However, both LaO and FeAs 
layers contribute approximately the same from -6 to -2 
eV, suggesting a hybridization between layers within this 
energy window. For AFM state, the calculations give a 
2.30 [is local magnetic moment on Fe by integrating the 
PDOS to Ep, which is similar to the value obtained from 
calculations for the bulk FeAs crystal. 
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TABLE I: Electron hopping tij and on-site energies t%j (in eV) 
matrix elements calculated from MLWFs. In AFM state, hop- 
ping will be different for iron atoms on different sub-lattices. 




FIG. 3: (Color online)Diagrams for (a) atomic-type MLWFs 
and (b) strongest hoppings in the system. The diagram is 
presented on the x-y plane with the most important hoppings 
labeled. Iron atoms and arsenic atoms are located at the 
vertices and centers of the square lattice, respectively. Note 
that (b) depicts an irreducible subset of hoppings, and the 
z-displacement of As atoms is not shown. 

To further understand the physics within the FeAs 
layers and connect our calculations to model calcula- 
tions, we have used the maximally localized Wannicr 
functions (MLWF) metho d 22 ! 23 to analyze the atomic 
orbitals which dominate the electronic structure near 
the Fermi surface. Sixteen MLWFs, including 10 d- 
type MLWFs on Fe and 6 p-type MLWFs on As, have 
been used to obtain a tight-binding effective Hamilto- 

nian H e ff = J2i e i c l c l + t i] c \ c ] by fitting the band 
structure around Ep. These MLWFs are then used to 
construct a model Hamiltonian matrix, from which we 
can regenerate the band structure using a tight-binding 
framework (blue curves in Fig. [T]) . In both AFM and PM 
cases, the tight-binding band structure fits the DFT band 
structure well, showing the validity of our model Hamil- 
tonian. We show the most important hopping terms in 
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Fig. [3l and the corresponding values are listed in ta- 
ble U together with the on-site energies. Due to the S4 
symmetry of the FeAs tetrahedra, the Fe 3d orbitals split 
into 3 non-degenerate (Z& x 2_ y 2 1 Sd^, 3d z 2) and 1 doubly 
degenerate energy state (3d a; (j / ) z ) in both PM and AFM 
states. Interestingly, the lowest lying 3d z 2 state in PM 
is the highest in AFM state, leaving the order of other 3 
states unchanged. The energy difference between lowest 
3d states and highest 4p orbitals (4p 2 ,(j / )) are ~ 1.2 eV 
and 1.6 eV in PM and AFM state, respectively, which 
is about the same magnitude as the Cu-0 splitting in 
cuprates. In the PM state, the strongest hoppings come 
from Fe 3da;2_j,2, 3d 2 2 and As 4p 2 orbitals, followed by 
the coupling from Fe 3d XZl 3d xy and As 4p^ orbitals. 
Remarkably, the direct hopping between As 4p x (y) i n the 
same x-y plane are as large as the Fe-As hopping. In 
the AFM state, since the mirror symmetry within the 
unit cell is removed, these hopping matrix elements split 
into two groups for spin-up electrons and spin-down elec- 
trons on Fe sites, respectively. Furthermore, the hopping 
between two neighboring As 4p z and 4p z orbitals are 
greatly enhanced in AFM states. The direct hoppings 
between Fe 3d orbitals are finite, but much smaller com- 
pared to Fe-As and As- As hopping in both cases. 

We finally present our calculations on doped LaOFeAs. 
The doping is simulated by substituting an oxygen atom 
in 8 primitive cells with a fluorine atom, and then fully 
relaxing the structure with the optimized lattice con- 
stants of the undoped system. Since each primitive 
cell contains 2 oxygen atoms, the doping corresponds to 
LaOi-aTzFeAs with x = 0.0625. The resulting material 
turns out to be still AFM, but with an extra total spin. 
With a;=0.125 doping, we reduced the number of prim- 
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FIG. 4: (Color online)DOS calculated from GGA+U with 
different Us, 2.0, 3.0, and 4.5 eV. E F is aligned to 0.0 eV in 
all cases. Due to degeneracy, only the a-spin DOS is plotted. 
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FIG. 5: (Color online)DOS of LaOi-^F^FeAs with a;=0.125 
projected onto FeAs layers. The Ef is aligned to 0.0 eV. 



AFM ground state has been found for undoped LaOFeAs 
via DFT calculations. We find that the geometry, elec- 
tronic structure and the magnetic state of this system are 
strongly related. In both AFM and PM states, the band 
structures around Fermi level are derived from Fe 3d and 
As 4p orbitals, and we have fitted bands crossing the 
Fermi surface to tight-binding Hamiltonians using ML- 
WFs. The parameters for the model Hamiltonians from 
the first-principles calculations can be used for modelling 
transport, magnetic and superconducting phenomena as- 
sociated with strongly correlated electrons in the system 
under investigation. While the system exhibits metallic 
behavior in DFT calculations, an inclusion of an on-site 
energy of 4.5 eV on Fe turns it into a semiconductor with 
a gap of 1.0 eV, which implies that the system is close to 
a Mott-type insulator. Due to the evident proximity of 
superconductivity to antiferromagnetism and the Mott 
transition, we suggest that the system may be a large- 
spin analog of the electron-doped cuprates, where AFM 
and superconductivity coexist. 

Note Added: Shortly after the first draft of this pa- 
per appeared, a linear spin density wave state was pre- 
dicted in an electronic structure calculation^, and dis- 
covered in neutron scattering experiments^, then exten- 
sively studied in [26|. We have compared the energy of 
such a magnetic state with the sublattice type AFM state 
discussed in this work, and can confirm that according 
to PWSCF calculations it is 109 meV per Fe lower in 
energy. 



itive cells involved by half, so that one oxygen atom in 
every 4 primitive cells was replaced with a fluorine atom. 
The fluorine substitution has an effect on electronic struc- 
ture which is primarily concentrated in the FeAs layers 
at energies near Ep, plus an impurity state at ~ 7.5 eV 
below Ep. We show the comparison of undoped and 
a;=0.125 doped system in Fig. The £=0.0625 doped 
system has a similar doping effect, but smaller in mag- 
nitude. The overall magnetic odering remains the same, 
but the magnetic moment is altered by 6%. 

In conclusion, we have performed first-principles cal- 
culations for LaOFeAs and LaOi-^F^FeAs systems. An 
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